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ABSTRACT 

We study the acceleration of the star HE 0437-5439 to hypervelocity and discuss its 
possible origin in the Large Magellanic Cloud (LMC) . The star has a radial velocity of 
723 kms -1 and is located at a distance of 61 kpc from the Sun. With a mass of about 
8 Mq, the travel time from the Galactic centre is of about 100 Myr, much longer than 
its main sequence lifetime. Given the relatively small distance to the LMC (18 kpc), 
we consider it likely that HE 0437-5439 originated in the cloud rather than in the 
Galactic centre, like the other hypervelocity stars. The minimum ejection velocity 
required to travel from the LMC to its current location within its lifetime is about 
500 km s -1 . Such a high velocity can only be obtained in a dynamical encounter with a 
massive black hole. We perform 3-body scattering simulations in which a stellar binary 
encounters a massive black hole and find that a black hole more massive than 10 3 M Q 
is necessary to explain the high velocity of HE 0437-5439. We look for possible parent 
clusters for HE 0437-5439 and find that NGC2100 and NGC 2004 are young enough 
to host stars coeval to HE 0437-5439 and dense enough to produce an intermediate 
mass black hole able to eject an 8 M Q star with hypervelocity. 

Key words: Stars: black-holes - Stars: individual (HE 0437-5439) - Methods: N-body 
simulations - Binaries: dynamics. 



1 INTRODUCTION 

Recent radial velocity measurements of young stars in the 
Galactic halo has led to the discovery of a new class of 
stars, the so-called hypervelocity stars (HVS), which travel 
with velocities exceeding the escape speed of the Galaxy 
jBrown et al.ll2005l , 120061 ). A total of 7 HVSs has been dis- 
covered so far, of which six are consistent with an ejection 
from the centre of the Milky Way. 

The existence of a population of HVSs was first pro- 
posed by Hills (1988), who considered hypervelocity ejec- 
tions as a natural consequence of galaxies hosting super- 
massive black holes (SMBHs). The interaction between a 
SMBH and a stellar binary can result in the dynamical 
capture of one of the binary components at the expense of 
the high velocity ejec t ion of its companion st ar (|Hillslll988l : 



lYu fc Tremaind 120031 : [Gualandris et alJl2005h . Later an al- 
ternative was proposed by Baumgardt et al.l (|2006l h who ar- 



gue that an intermediate-mass black hole (IMBH) can eject 
stars with extremely high velocities as it spirals in towards 
the Galactic centre. The latter mechani sm tends to pro- 
duce HVSs in bursts which last a few Myr l|Baumgardt et all 
2006). Since the reconstructed ejection times of the observed 
HVSs are roughly evenly distributed over about 100 Myr 
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l|Brown et al.|[20 0r?). it seems likely that they were ejected 
in a continuous process. For 6 of the observed HVSs, in fact, 
the travel time from the Galactic centre and the estimated 
age are consistent and, by assuming a suitable proper mo- 
tion, it is possible to calculate a realistic trajectory. 

For HE 0437-5439, however, an ejection from the Galac- 
tic centre appears problematic as the time required to travel 
from the centre to its observed location exceeds the maxi- 
mum lifetime of the star. HE 0437-5439 is a main-sequence 
star with spectral type B and a mass of about 8.4 Mq for a 
solar metallicity or of about 8 Mq for a, Z = 0.008 metallic- 
ity, which corresponds to main-sequence lifetimes of 25 and 
35 M yr, respectively. Adopt ing a distance of 61 kpc from the 
Sun IjEdelmann et al] 120051 ) and considering a total space 
velocity equal to the radial velocity only, the travel time for 
a straight orbit from the Galactic centre is of about 80 Myr, 
much longer than the main-sequence lifetime of the star. 
The estimated proper motion for which the orbit of the star 
comes within 10 pc from the centre is of about 0.55 masyr^ 1 
jEdelmann et al.ll2005l h This results in a tangential velocity 
of about 160kms _1 and hence in a total space velocity of 
740 km s" 1 , which corresponds to a travel time of 100 Myr 
for a realistic orbit. 

We discard the possibility that HE 0437-5439 is a blue 
straggler originated in the Galactic centre since, even in the 
case of a merger between two less massive stars, the merger 
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product would not survive as an 8 Mq star for much longer 
than the main-sequence lifetime of such a star. 

Given the fact that HE 0437-5439 is much closer to 
the L MC galaxy than to the Milky Way, lEdelmann et al.l 
( 2005) suggest that HE 0437-5439 might have been ejected 
from the LMC. Considering a total distance of 18kpc from 
the LMC centre, the minimum ejection velocity required to 
travel to the current position is about 500kms _1 . (Here 
we adopted a travel time of 35Myr, which is equal to 
the main-sequence lifetime in a low metallicity environ- 
ment.) The radial velocity of HE 0437-5439 relative to the 
LMC is 461 km s -1 , and the required tangential velocity is 
then 160kms _1 , which in turns implies a proper motion of 
1.9masyr _1 . A velocity of the order of 500kms _1 can only 
be obtai ned in a dynamica l intera ction with a massive black 
hole (see iGualandris et al.l (I20051 ) 1 ). 

If HE 0437-5439 was ejected by an IMBH and if it came 
from the LMC, the natural consequence is that there must 
then be an IMBH in the LMC. We investigate this hypoth- 
esis by studying the dynamical ejection of an 8 Mq main- 
sequence star in an interaction with a massive black hole. 
We simulate the encounter between a stellar binary and a 
black hole to study the minimum mass of the black hole 
required to accelerate one of the interacting stars to a veloc- 
ity of at least 500kms _1 . We then investigate where in the 
LMC such an interaction could have taken place. 

We find that a black hole of > 10 3 M is re- 
quired to explain the velocity of HE 0437-5439. Such an 
IMBH could form in a y oung and dense star cluster 
l|Portegies Zwart et al ]|2004 ) and would help driving the fre- 
quent strong encounters with massive stars needed to make 
high- velocity ejections likely (Baumga rdt et all [20061. The 
IMBH must still be present in a star cluster which con- 
tains stars as massive as HE 0437-5439. In addition, the 
parent star cluster must have been massive enough and 
with a short enough relaxat ion time to produce an IMBH 
l|Portegies Zwart et al.l [20041 ). The most likely clusters are 
NGC 2100 and NGC 2004. 



2 THREE-BODY SCATTERINGS WITH A 
MASSIVE BLACK HOLE IN THE LMC 

In this section we explore the hypothesis of a dynamical 
ejection from the LMC by means of numerical simulations of 
three-body scatterings with a massive black hole. In § !2.1l we 
focus on interactions in which a binary containing a young 
8Mq main-sequence star (representing HE 0437-5439) en- 
counters a single black hole of 10 2 Mq to 10 4 Mq , whereas 
in § 12.21 a single 8 Mq main-sequence star encounters a bi- 
nary black hole. In the first case, an exchange interaction 
can lead to the high velocity ejection of one of the binary 
components while in the latter case the main-sequence star 
can be accelerated in a fly-by. 

The simulations are carried out using the sigma3 pack- 
age, which is part of the STARLAeQ] software envir onment 
([McMillan fc Hutll 19961 ; iPortegies Zwart et al.ll200ll ). 

For each simulation we select the masses of the three 
stars, the semi-major axis and eccentricity of the binary and 
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the relative velocity at infinity between the binary's centre 
of mass and the single star. The phase of the encounter- 
ing binary is random ly drawn from a uniform distribution 
jHut fc Bahcalllll983l ) and the orbital eccentricity is taken 
randomly between circular and hyperbolic from the thermal 
distribution. The impact parameter b is randomised accord- 
ing to b 2 in the range [0 — 6 ma x] to guarantee that the prob- 
ability distribution is homogeneously sampled. The maxi- 
mum value femax is determined automatical ly for each set 
of experiments (see IGualandris et al"] |2004) for a descrip- 
tion). Energy conservation is typically better than one part 
m 10 6 and, in case the error exceeds 10 5 , the encounter is 
rejected. The accuracy in the integrator is chosen in such a 
way that at most 5% of the encounters are rejected. Dur- 
ing the simulations we allow for physical collisions when the 
distance between any two stars is smaller than the sum of 
their radii. The black hole is assumed to be a point mass, 
while for t he stars we adopt zer o-age main-sequence radii, 
taken from lEggleton et all l|l989l ). 

2.1 Encounters between a stellar binary and a 
massive black hole 

We consider encounters between a binary consisting of two 
main-sequence stars with masses mi and mi and a single 
black hole with mass Mbh- The mass of one binary compo- 
nent is fixed to 8 Mq while the mass of the companion star 
has values of 2 Mq , 4 Mq , 8 Mq and 16 Mq . For each choice 
of stellar masses, the semi-major axis a is varied between 
a minimum value, which is set by the physical radii of the 
stars so that the two components do not touch at the first 
pericentre passage, and a maximum of 1AU. We consider 
black holes of masses 10 2 M , 10 3 and 10 4 M . The relative 
velocity at infinity between the black hole and the binary's 
centre of mass is set equal to 20kms~ 1 . 

In Fig.Q] we present the probability of different out- 
comes (branching ratios) in the simulations of an encounter 
between an equal mass binary (mi = m,2 = 8M ) and a 
10 3 M black hole. For each of the initial semi-major axes 
we perform a total of 1500 scattering experiments, which re- 
sult either in a fly-by, an exchange or a merger. Mergers oc- 
cur preferentially in the case of tight binaries, leaving place 
to preservations and exchanges for wider binaries. In an ex- 
change interaction one of the main-sequence stars is cap- 
tured by the black hole while the other star is ejected, pos- 
sibly with high velocity. If the velocity exceeds 500 km s~ , 
we regard the star as a possible candidate for HE 0437-5439. 
This occurs in about 10% or less of all encounters for semi- 
major axes in the range 0.1 — 1.0 AU. 

In Fig. [2] we show the branching ratios for high velocity 
ejections for binaries with a 8Mq star and a companion 
with mass in the range 2.0— 16.0 Mq. While the 8Mq main- 
sequence star is ejected, the companion star is captured by 
the IMBH. As in the previous plot, the black hole mass is 
of 10 3 Mq . Exchange encounters with fast escapers are more 
likely in the case of binaries with larger companion masses 
and/or shorter orbital periods. 

Analytical estimates by lYu fc Tremainel l|2003h for the 
case of tidal breakup predict an ejection velocity at infinity 
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Figure 1. Branching ratio for the outcome of encounters between 
a stellar binary and a single 10 3 M black hole as a function of 
the initial binary semi-major axis. The two binary components 
are assumed to be 8Mq main-sequence stars. The different out- 
comes are: merger (stars), preservation (triangles) and exchange 
(squares). The subset of exchanges with a high velocity escaper 
(Voo ^ 500 kms -1 ) is indicated with bullets. The error bars repre- 
sent the formal (ltr) Poissonian uncertainty of the measurement. 
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Figure 2. Fraction of exchange encounters between a stellar bi- 
nary and a 10 3 Mq black hole that results in the ejection of a 
main-sequence stars with a velocity ^ 500 kms -1 . The bullets 
indicate the relative probability for an incoming binary with two 
8M0 stars (see Fig.[T]l. The other curves are computed with dif- 
ferent companion masses to an 8Mq star: 2Mq (crosses), 4Mq 
(squares), 16 Mq (triangles). In the latter case, the ejected star 
is the secondary. 



where v is a function of the distance of closest approach 
.Rmin- If we define a dimensionless closest approach param- 
eter 7?„;„ = flmi " 1,/3 , v varies in the ranee 130- 

mm a \ m / ' ° 

160 kms -1 for = 0-1. For m = 8 Mq and a = 0.1 AU 
Eq-EDyields v,^ w 560 kms" 1 for M hh = 10 2 M Q ; a 100 M 
IMBH would in principle be sufficient to explain the origin 
of HE 0437-5439. 

In Fig. [3] we show the velocity distributions of the es- 
caping single star in the case of equal mass binaries (mi = 
rri2 — 8 Me). The different panels refer to different values of 
the black hole mass, from 10 2 M Q (left) to 10 3 M (middle) 
and 10 4 Mq (right). In each panel, the different distribu- 
tions refer to three different values of the initial semi-major 
axis. While it appears possible to reach the required recoil 
velocity in the case of Mhh = 10 2 Mq , the smallest possible 
semi-major axis is needed (a = 0.1 AU) for this to happen, 
and only in about 10% of exchange encounters, which corre- 
sponds to about 3% of all simulated scatterings. In the case 
of the 10 4 Mq black hole it is possible to achieve such high 
velocities even for rather wide binaries. Based on these data 
we conclude that a 100 Mq black hole is unlikely to have 
resulted in the ejection of HE 0437-5439 and a black hole 
mass £ 10 3 Mq is favoured for typical values of a. 

We note that the velocity distributions are not sensitive 
to the initial velocity between the binary and the single star 
as the total energy of the system is dominated by the bind- 
ing energy of the binary rather than by the kinetic energy of 
the incoming star. Additional experiments performed with 
initial velocities of 5 kms -1 and 10 kms -1 showed no appre- 
ciable difference in the ejection velocities. 

We also note that for detached binaries the binary tidal 
radius is always larger than the tidal radius of the single stars 
and therefore tidal breakup and ejection occur before the 
stars can be disrupted by the black hole. For (near) contact 
binaries, instead, it is possible to disrupt a star before an 
ejection can take place. This process does not depend on the 
black hole mass but on the size of the binary compared to 
the size of the stars. Our simulations only include detached 
binaries as the proper treatment of contact binaries would 
require hydrodynamical simulations. 

The analytical estimates obtained with Eq.[T]are consid- 
erably larger than the results of our simulations (see Fig.[3j) . 
This discrepanc y is pr obably caused by the assumption of 
lYu fc Tremainel (|2003l ) that the binary is disrupted well 
within the tidal radius, i.e. R min ^ 1, whereas in 10% to 
20% of our simulated encounters i? min > 1. For example, for 
R min — 3 — 4, a capture can occur which eventually leads to 
the ejection of one star. Since the ejection velocity depends 
largely on the Keplerian velocity of the binary at the posi- 
tion of tidal breakup, we expect a dependence of the final 
velocity of escapers on the parameter J? m i n . Fig. [4] shows the 
velocity at infinity of the ejected star versus the dimension- 
less distance of closest approach. In this example we consider 
equal mass binaries with m — 8 Mq and a black hole of mass 
1O 3 M . The fi gure shows that the highest ejection veloci- 
ties are preferentially achieved in cl ose encounters. The ne - 
glect of relatively wide encounters bv lYu fc Tremainel (120031 ) 
might therefore be the origin of the apparent discrepancy in 
the velocities. The agreement betwee n our numerical simula - 
tions and the analytical estimates of lYu fc Tremainel (|2003l ) 
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Figure 3. Velocity distributions for different values of the black hole mass M hh = 10 2 M Q (left), 1O 3 M (middle), 10 4 M Q (right), and 
of the initial semi-major axis a = 0.1, 0.3, 1.0AU. 



1500 




V_ (km s" 1 ) 



Figure 4. Velocity of the ejected star after an exchange inter- 
action with a 10 3 Mq black hole as a function of the distance of 
closest approach to the black hole. The different symbols indicate 
different initial values of the binary semi-major axis. 



Figure 5. Velocity distribution for the ejected star after an in- 
teraction between an equal mass (m = 8Mq) binary and a black 
hole of 1O 2 M (black), 10 3 M (blue)and 10 4 M (red). These 
velocity distributions are integrated over the entire range of or- 
bital separations for the initial binary. 



improves for higher black hole masses, and is very go od for 
encounters with a SMBH (see I Gualandris et al.l r2005). 

In our systematic study of the effect of the initial semi- 
major axis of the interacting binary we adopted a homoge- 
neous sampling in log a. Furthermore, the number of scat- 
tering experiments performed per initial selection of a are 
weighted with equal cross section. If the distribution of or- 
bital separations in a star cl uster is flat in log a, like in the 
case of young star clusters (jKouwenhoven et al.l 1200*51 ) . we 
can superpose the results of these experiments in order to 
acquire a total velocity distribution of the ejected star. In 
Fig. [5] we present this superposed velocity distribution for a 
binary consisting of two 8 M stars. The three histograms in 
this figure give the velocity distribution for a 100 M (left), 
10 3 M (middle) and 10 4 M (right) black hole. 



2.2 Encounters between a single star and a binary 
black hole 

We now consider encounters betwee n a single star and a 
binary black hole. In a recent study iGiirkan et alj (|2006t ) 
argue that it is possible to form a black hole binary as a 
result of two collision runaways in a star cluster. The bi- 
nary eventually merges due to the emission of gravitational 
wave radiation (|Peterslll964 ) but, before coalescing, the two 
black holes experience frequent interactions with other clus- 
ter stars, possibly accelerating them to larger velocities. 

We simulate these encounters and study the probability 
of high velocity ejections as a function of the binary param- 
eters: the masses of the two black holes Mbm and Msh2 
and the initial semi-major axis a. The mass of the single 
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star is fixed to 8M , in accordance with the observed mass 
of HE 0437-5439. We consider black hole masses of 10 2 M , 
10 3 M and 10 4 M in both equal mass and unequal mass 
binaries. For each choice of black hole masses, the semi- 
major axis is varied between a minimum value of 1 AU and 
a maximum value of 1000 AU. The minimum value is cho- 
sen in order to guarantee that the binary remains unaffected 
by the emission of gravitational wave radiation during the 
period over which the encounter takes placfl 

Except for a few mergers (about 1-2% of all cases), the 
vast majority of the encounters result in a fly-by. The ve- 
locity at infinity of the escaping star depends sensitively on 
the impact parameter of the encounter. Only in encounters 
for which the impact parameter is comparable to the binary 
separation are high velocity ejections realized. Most of these 
encounters are resonant, allowing the incoming star to come 
as close to one the black holes as a fraction of the orbital 
separation. Only a few % of the encounters result in ejection 
velocities as high as 500 kms -1 and we argue that this mech- 
anism is unlikely to have produced the observed velocity of 
HE 0437-5439. 



3 AN IMBH IN THE LMC? 

The analysis carried out in the previous sections suggest the 
intriguing idea that HE 0437-5439 was ejected by an en- 
counter with an IMBH of > 10 3 M in the LMC. Such 
an IMBH would most likely be found in a young dense 
cluster containing s tars c oeval to HE 0437-5439. Recently 
iMackev fc Gilmorel (|2003l ) measured the structural param- 
eters for 53 LMC clusters. A total of nine clusters in this 
database are younger than 35 Myr. Three of these clusters, 
NGC 1818, NGC 1847 and NGC 1850, have the appropriate 
age to be the host of HE 0437-5439 but their current den- 
sities are much smaller than what is required to produce 
an IM BH |Portegies Zwart et all 12004 iMiller fc Hamilton! 
12002! '). Of the six remaining clusters, three (NGC 1805, 
NGC 1984, NGC 2011) have a mass below 5000 M Q and are 
therefore unlikely to have produced an IMBH. Also, the 
presence of an IMBH in these cluste rs would have been no- 
ticed easily dBaumgardt et al1l2005h . 

The remaining three clusters are listed in Tab.[T] For 
these clusters we recons truct the initial relaxation time via 
the method described bv lPortegies Zwart fc C hen (200^). In 
this relatively simple model, the current relaxation time has 
a one-parameter relation with the cluster's initial relaxation 
time. Of the three clusters listed in Tab.[TJ one (R 136) is 
too young to produce an IMBH and experience a strong en- 
counter with a binary. The best candidate clusters to host 
an IMBH are then NGC 2100 and NGC 2004. We can es- 
timate the mass of a possible black hole in these clusters 
by adopting the relation between the cluster structural pa- 
rameters and the m ass of a central black hole presented by 
iHeggie et al.l {2006). The authors performed direct N-body 
simulations to quantify the relation between the mass of a 
central IMBH and the ratio of the core radius to the half- 
mass radius. The last column in Tab. [T] provides an estimate 



Table 1. List of young (age ,$ 35 Myr) LMC clusters with 
an estimated initial relaxation time smaller than 100 Myr. The 
columns give the name of the cluster, the mass, the core radius, 
the ratio of the core radius to the half-mass radius, the age, the 
current two-body relaxation tim e, the initial two-body rel a xation 
time (computed using Eq. 1 of iPortegies Zwart fc Chenl d2006l) 
with k = 0.1) and the black hole mass based on lHe^gieet**aU 
(2006). For estimating the relaxation time we adopted a mean 
mass (m) = 1 and a Coulomb logarithm parameter 7 = 0.4 for 
the lower limit, and (m) = 0.65 and 7 = 0.01 for th e upper limit. 
The c ore and half-mass radius for R 136 are from iBrandl et al.l 
lll996h . 



name 


M 

M 


r c 
pc 


r c /r h 


age 
Myr 


Myr 


T o 

rlx 

Myr 


M 


R136 


35000 


0.32 


0.29 


3 


7-17 


6-14 


1000 


NGC 2004 


27000 


1.57 


0.50 


20 


71-170 


41-98 


1600 


NGC 2100 


30200 


1.22 


0.62 


16 


51-120 


31-73 


2200 



for the mass of the IMBH using this relation. We find that, 
given the structure of the observed clusters, a 1600 M to 
2200 M IMBH could be present. 

We argue that NGC 2100 and NGC 2004 form the most 
promising birth places for HE 0437-5439. If this were the 
case, the travel time from the parent cluster to the cur- 
rent location would be less than about 20 Myr. Such a short 
travel time would require an ejection velocity ^ 800 kms 
for the adopted total distance to the LMC of 18kpc. This 
large velocity could only be obtained with a black hole of 
several 10 3 M , somewhat higher than the estimated mass 
of a possible black hole in these clusters. 

We now compute the rate for ejection of hyperveloc- 
ity stars from NGC 2100 and NGC 2004 using the param- 
eters listed in Tab.[T] and a dimensionless cross-section for 
exchange encounters a ~ 5000, as derived from the scatter- 
ing experiments. Adopting a central density of 2 x 10 4 pc -3 
(consistent with a King Wo = 9 initial profile) , a velocity of 
10 kms" 1 and a semi-major axis of 0.2 AU, we obtain an en- 
counter rate of about one per 2 Myr for both clusters. Since 
only one in about ten exchange encounters produces a fast 
escaper (see Fig. [2]), we derive a production rate of one per 
20 Myr. This value represents a lower limit as we don't take 
into account the effects of a mass function in the core and 
the semi-major distribution. It is therefore conceivable that 
NGC 2100 or NGC 2004 has produced one high velocity es- 
caper, in which case an IMBH must be present in one of 
these clusters. 
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2 Two 10 3 M black holes in a 1 AU orbit will merge due to the 
emission of gravitational waves in about 24 Myr. 
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